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'  The  overall  stfective  thermal  resistance  of  a  corrugated 
fiberboard  container  was  determined  as  a  function  of 
air  velocity  and  board  thickness,  end  an  appropriate 
design  curve  is  presented.  The  thermal  resistance  of 
the  container  is  treated  as  a  sum  of  resistsnces,  with 
individual  resistance  values  presented  tor  the  interior 
interfacial  contact  resistance,  the  board  resistance, 
and  the  exterior  boundary  resistance  Behavioral  dif 
ferences  were  found  for  heating  versus  cooling  at  con¬ 
ventional  board  thicknesses  and  a  possible  explanation 
for  this  behavior  is  presented 


United  States 

WPIf«nPn  OT 

Agriculture 


Fond 

Products 

Laboratory' 

Mseaerch 

F??*40e 

August  itot 


Overall  Effective 
Thermal  Resistance 
of  Corrugated 
Flberboard 
Containers 


•y 

OAVIO  W  BOMMETT,  Chemical  Engineer 


Ovsr  tits  yoars  ol  its  development,  the  properties  of  cor- 
rugeted  flberboard  have  been  studied  extensively  in 
order  to  predict  Its  behavior  in  usage.  As  ihees  proper¬ 
ties  became  better  understood  the  corrugated  container 
took  the  place  of  the  wood  boa  for  general  use  Packag¬ 
ing  has  now  evolved  to  the  point  where  the  corrugated 
container  is  the  primary  packaging  material 

However,  despite  the  generation  of  a  broad  data  base 
on  the  mechanical  properties  of  corrugated  flberboard, 
little  indapth  research  has  been  undertaken  on  its 
nonmechanical  properties  m  particular,  even  though 
the  agricultural  induet ry  ships  vast  quantities  of  pro¬ 
duce  in  oorrugated  boxes,  little  knowledge  is  available 
on  the  thermal  properties  of  the  oorrugated  container. 
White  shipping  and  storing  produoe  and  frocen  foods 
has  been  studied  extensively,  research  emphasis  has 
generally  bean  on  the  carried  product,  with  only  mild 
consideration  given  to  the  package  itself.  Consequent- 
ly.  it  is  difficult  to  extract  from  the  literature  any  design 
criteria  tor  the  thermal  properties  of  corrugated  boxes. 

Hamaher*  dM  determine  the  thermal  reeletanoe  and 
conductivity  of  ovendry  board  as  a  function  of  board 
thickness  and  components  and,  baaed  on  these  results, 
he  generated  a  mathematical  mods!  to  prediet  its  ther¬ 
mal  reaiatanoa.  While  the  thermal  reeletanoe  values  are 
oorrect,  the  use  of  these  values  alone  proved  inade¬ 
quate  tor  predicting  the  actual  thermal  performance  of 
bases  In  see.  Other  controlling  factors  weeded  to  be 
coupled  wMh  the  thermal  reeietence  to  describe  the  box 
behavior. 


As  e  result,  this  study  wee  undertaken  to  determine  the 
overall  effective  thermal  reeletanoe  of  a  oorrugated 


flberboard  container  The  factors  considered  were 
board  thickness,  air  velocity  around  the  box.  and  direc 
tton  of  the  thermal  flux  (inward  or  outward  heat  flow 
per  unit  area  per  unit  time) 


Techniques  exist  for  predicting  the  time-temperature 
behavior  of  s  product  packed  in  a  corrugated  flberboard 
box.  The  user  of  these  techniques  is  required  to  know 
(t)  the  geometry  of  the  product.  (2)  its  thermal  charac 
t eristics  of  specific  heat,  density,  and  thermal  conduc¬ 
tivity;  (3)  the  external  bulk  air  environment  in  which  the 
heating  or  cooling  occurs,  and  (4)  the  thermal  behavior 
between  the  product  end  the  bulk  air  it  is  first 
necessary  to  understand  the  conditions  satisfying  the 
first  three  requirements  in  actual  use  and  to  see  how 
these  conditions  were  simulated  experimentally 

Wsgulremsato  U1 

While  the  thermal  and  geometric  properties  of  s  pro¬ 
duct  are  critical  to  any  analysis  (requirements  i  and  2). 
the  diversity  of  products  is  so  extensive  that  any 
amount  of  consideration  of  these  variables  would  still 
be  inoompiete.  Therefore,  in  this  study  product 
variables  were  kepi  constant  by  smptoytng  only  one 
commercial  product  in  one  standard  sue.  namely 
464-gram  flat  packs  of  margarine. 

For  requirement  9,  from  the  user's  standpoint  two  com¬ 
mon  bulk  air  environments  are  generally  encountered 
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(1)  Th«  package  at  a  given  temperature  is  placed  in  a 
constant-temperature  room  at  a  different  temperature 
and  at  a  very  low  circulating  air  flow  rate  such  that  the 
package  is  subjected  primarily  to  convection  currents; 

(2)  high-velocity  air  at  a  temperature  different  than  the 
package  strikes  the  package  on  a  single  face  and  flows 
around  it.  A  third  situation  frequently  arising  is  air  flow¬ 
ing  through  boxes  by  means  of  vent  holes.  This  special 
case  is  being  evaluated  as  a  separate  problem  to  be 
reported  on  at  a  later  date.  In  this  study  the  air  flow 
patterns  used  were  intended  to  simulate  the  above  two 
common  environments. 


Having  set  the  conditions  for  the  first  three  re¬ 
quirements.  one  can  then  consider  the  fourth  require¬ 
ment  A  commercial  packer  or  user  is  rarely  concerned 
with  the  container  but  only  with  the  product  inside.  In 
this  respect  the  container  serves  as  part  of  the  overall 
resistance  to  the  flow  of  heat  between  the  product  and 
bulk  air 

Consider  the  package  as  being  three  separate  thermal 
resistance  zones  through  which  heat  must  pass:  (1)  An 
mterfacial  contact  resistance  between  the  surface  of 
the  product  and  the  inner  surface  of  the  corrugated 
ftberboard;  (2)  the  resistance  of  the  corrugated  fiber- 
board  itself;  and  (3)  the  resistance  between  the  outer 
surface  of  the  ftberboard  and  the  bulk  air.  The  sum  of 
these  three  resistances  is  the  overall  effective  thermal 
resistance  of  the  package  Its  reciprocal  is  the  overall 
effective  heat  transfer  coefficient 

To  satisfy  the  fourth  requirement  a  user  must  know  the 
overall  resistance,  which,  unfortunately,  changes  with 
both  time  and  point-by-point  location  on  the  box.  The 
changes  are  due  to  the  heterogeneous  contact  between 
box  and  product,  with  the  thickness  of  the  ftberboard 
on  each  side  of  the  box.  and  with  the  turbulent  air  pro¬ 
file  surrounding  the  box. 

While  a  strict  analysis  would  take  each  of  these  factors 
into  account,  the  results  would  be  overly  cumbersome 
for  general  field  use.  Thus,  this  work  is  an  attempt  to 
simplify  the  overall  time-position  resistance  history  into 
a  single-value,  approximate  resistance  characterizing 
the  overall  effective  thermal  behavior  of  the  package. 
For  additional  theoretical  considerations  see  Appendix  I. 

The  actual  technique  for  using  the  resultant  resistance 
values  for  field  use  prediction  of  a  time-temperature 
history  on  an  arbitrarily  sized  container  having  an  ar¬ 
bitrary  product  will  be  published  in  a  subsequent 
report. 


tats  of  four  464-gram  (g)  packs  of  commercial  flat  pack 
margarine  were  stacked  in  a  130-  x  130-  x  130-millimeter 
(mm)  cube  and  packaged  in  corrugated  ftberboard  con¬ 
tainers  Board  thickness  ranged  from  0.33  mm  (liner- 
beard  only)  to  81  mm  of  built-up  corrugated  ftberboard 
The  various  packages  were  subjected  to  heatlng/ooot- 


Tsble  1.— Overall  effective  heat  transfer  coefficients  and  ther¬ 
mal  resistances  of  a  cubical  corrugated  container 
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Heat 

transfer 
coefficients  (U) 
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0 
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0.13 
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ing  periods  at  bulk  air  velocities  of  0,  1.8,  3.6.  and  8.4 
meters  per  second  (m/s).  Instantaneous  overall  heat 
transfer  coefficients  were  determined  at  the  center 
point  of  each  side  of  the  package  at  frequent  intervale 
for  the  duration  of  the  temperatury  change.  For  each 
environmental  condition,  an  overall  time  and  position 
average  effective  heat  transfer  coefficient  (U)  was 
calculated  for  each  package  as  a  whole.  The  resulting 
coefficients  are  presented  in  table  1  along  with  their 
reciprocals  (j).  which  are  the  desired  overall  effective 
thermal  resistances  (R)  (fig.  1).  For  the  detailed  Ex¬ 
perimental  Procedure,  see  Appendix  II. 


Satisfying  requirement  3,  the  data  (fig.  1)  for  no  air  flow 
(maximum  reeistanca,  line  III)  and  Infinite  air  flow 
(minimum  resistance,  line  II)  represent  the  primary 
design  data.  Values  tor  1.8,  3.6,  and  6.4  m/s  are  omitted 
as  they  are  sufficiently  cfose  to  the  Infinite  sir  velocity 
values. 

The  throe  llnee  of  figure  1  represent  three  thermal 
reeistanoe  levels:  The  ovendry  board-only  reetstanoe 
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Table  1.— Overall  effective  heat  transfer  coefficients  and  ther- 
mal  resistances  el  a  cubical  corrugated 
container— oon. 


Air 

velocity 

Board 

thickness 

Heat 

transfer 
coefficients  (U) 

Thermal 

resistance 

(R  -  11V) 

m/9 

mm 

W/K»m* 

K-m’/W 

HEATINO 

0 

0 

3.06 

0.33 

33 

2.38 

.42 

423 

2.21 

.45 

12.7 

1.64 

.61 

254 

1.25 

.80 

50.8 

.91 

1.10 

1.8 

0 

7.21 

.14 

33 

5.33 

19 

4.23 

3.23 

.31 

12.7 

1.99 

.50 

254 

1.36 

.73 

508 

.96 

1.04 

36 

0 

670 

.15 

33 

5.33 

.19 

4.23 

3.46 

29 

12  7 

2.04 

.49 

25.4 

1  36 

.73 

50.8 

.96 

1.04 

9.4 

0 

783 

.13 

33 

5.50 

18 

423 

3.63 

28 

12  7 

190 

.50 

25.4 

1.42 

.70 

508 

96 

1.04 

• 

0 
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.13 

33 

5.56 

18 
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28 
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.50 
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1.36 
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50.8 

98 

1.04 
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(line  I),  the  board  plus  board-product  interfaclal 
resistance  (line  II),  and  the  total  overall  resistance 
which  includes  the  maximum  outer  bulk  air-board  inter- 
facial  resistance  (line  III). 

Line  I  is  taken  from  the  equation1  for  ovendry  cor¬ 
rugated  fiberboard  at  a  mean  temperature  of  -  1  *  C 
corresponding  to  the  experimental  mean  temperature  of 
this  study.  Lower  resistance  values  than  shown  in  line  I 
may  be  in  order  due  to  moisture  (5  to  10  pet)  present  in 
the  board  during  the  experimental  runs.  The  extent  of 
the  deviation  is  not  accurately  known,  but  Ramaker 
presented  evidence  of  a  possible  15  to  80  percent 
decrease  in  resistance  due  to  moisture.  This  effect  is 
currently  being  investigated. 

Line  II  represents  the  extrapolated  zero-intercept  values 
from  table  2.  These  values  are  for  infinite  bulk  air 
velocity,  implying  no  outer  bulk  air-board  boundary 
resistance  exists.  Only  the  board  resistance  and  the 
board-product  interfacial  resistance  remain.  The  shape 
of  line  II  should  be  noted.  Because  the  box  and  the 
solid  have  a  tight  fit,  one  would  not  expect  convection 
at  the  interface.  Nor  should  radiation  be  a  significant 
factor  due  to  the  small  temperature  difference  across 
the  interface.  Thus  it  would  seem  reasonable  to  expect 
a  constant  contact  resistance  at  this  interior  interface, 
demonstrated  by  the  parallelism  of  lines  I  and  II.  This  is 
not  the  case  as  is  emphasized  by  figure  2  (line  II  minus 
line  I)  which  represents  the  interior  interfacial  contact 
resistance  only.  II  remains  to  be  determined  if  there  is 
a  true  unexplained  mechanism  taking  place  or  if  the  ef¬ 
fect  Is  merely  scatter  in  the  data.  Regardless,  one  can 
assume  that  the  board-product  interfacial  resistance 
for  this  system  will  be  0.18  to  0.27*  Kelvin  •  square 
meters  per  watt  (K^m'/W)  at  normal  box  thicknesses1  of 
2.5-20  mm.  This  implies  a  minimum  overall  combined 
resistance  of  three  times  that  due  to  board  alone.  This 
means  that  the  board  resistance  alone  is  not  sufficient 
for  calculating  heating/cooling  times  for  boxes  and  the 
inner  and  outer  interfacial  resistances  must  be  in¬ 
cluded. 

Figure  1,  line  III  represents  the  overall  effective 
resistance  for  no  bulk  air  flow  (stagnant  air),  that  is,  ex¬ 
terior  convection  effects  only.  The  outer  air  boundary 
layer  would  have  minimal  disturbance,  therefore  max¬ 
imum  resistance.  As  with  the  difference  between  lines  I 
and  II,  the  difference  between  lines  II  and  III  is  the  ac¬ 
tual  thermal  resistance  of  the  outer  boundary  air  layer. 
For  board  thicknesses  greater  than  20  mm  the  stagnant 
outer  air  resistance  is  nearly  constant  at  0.00  K«mVW. 
However,  in  the  normally  used  thickness  range  of  less 
than  20  mm  an  obviously  differing  pattern  exists. 
Heating  versus  cooling  shows  a  distinct  difference  as 
shown  by  the  splitting  of  line  III.  One  possible  explana¬ 
tion  for  this  behavior  can  be  found  by  considering  the 
nature  of  the  outer  boundary  air  layer.  A  hypothesis  of 
this  behavior  can  be  found  In  Appendix  III. 


1  Typical  UMcfcnass  rang**  fof  single-wall  combined  board  are  4  6-58 
mm  (A-ftuto,  3.S-4.I  mm  (C  flutci,  and  28-31  mm  (B-tlute). 
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Table  2.— Determination  of  effective  thermal  reeietanoe  tram  air 
«atoetty'J 


Flux 

Raeiatanca  (K«m,/W) 

•  a  *  tbW 

direction 

a 

b 

mm 

0.33 

AaaliaaAaalUa 

c 

0.21 

0.0023 

(kraft 

liner) 

h 

.18 

.0206 

4.2 

c 

.27 

.0026 

(single 

wall) 

h 

.28 

.0144 

127 

c 

.50 

.0049 

h 

.50 

.0096 

25.4 

c 

.77 

.0068 

h 

.73 

.0085 

50.8 

c 

1.10 

.0006 

h 

1.04 

.0027 

'  Constant  "b"  ganarated  from  linear  regression  of  data  of 
table  1.  excluding  data  for  infinite  air  velocity.  Extrapolation 
of  theae  regression  lines  to  zero  inverse  air  velocity  gives  the 
constant  "a,"  corresponding  to  the  infinite  air  velocity  data  in¬ 
cluded  in  table  1. 

’  Equations  are  valid  for  air  velocities  in  the  range  0.09  to  5.4 
m/s.  For  velocities  in  the  range  0.00  to  0.09  m/s  use  data  for  0 
m/s  in  table  1  or  figure  1. 

1  The  value  "x"  is  the  air  velocity  in  meters  per  second. 
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While  the  proposed  mechanism  needs  to  be  proved  or 
disproved,  the  varying  therms!  behavior  does  exist  In 
the  region  of  common  board  thicknesses.  Most 
significantly,  the  thermal  resistance  for  a  single-wall 
box  during  heating  with  no  air  movement  except 
natural  convection  is  0.44  K«m*/W  or  five  times  the 
resistance  of  board  alone.  Again,  this  emphasises  the 
necessity  of  Including  the  Interfacial  resistances  in  any 
calculations.  Similarly  with  no  forced  air  movement  a 
thawing  (heating)  period  may  take  considerably,  longer 
than  a  freezing  (cooling)  period  as  shown  by  the  higher 
reeistance  (line  lll-a)  for  heating  compared  to  the  lower 
reelstance  (line  lll-b)  for  cooling. 
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Figur *  3  —  C hong*  with  tim*  of  th*  overall  effective  h**t  tr*n»t*r  co*tt> 
c root  (U>  *$  reflected  by  th*  tint  root  0}  of  th a  equation 
U/h  *  k  fn  i L  I A}  cooling.  (&}  hooting  MoriiOnin  tin*  it  th* 
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At  normal  board  thicknesses  blast  cooling  of  a  box 
with  no  vent  holes  will  have  only  marginal  effects.  Com 
pare  lines  lll-b  and  II  In  figure  1.  Conversely,  blast 
heating  would  be  a  definite  asset  to  reducing  heating 
or  thawing  times  (line  lll-a  minus  line  II). 

Noting  that  the  values  of  resistance  reported  here  are 
average  values  calculated  for  temperature  changes  90 
percent  or  more  complete,  examination  of  figure  3 
shows  that  If  one  were  to  abbreviate  the  heating/cool¬ 
ing  period,  the  average  value  of  A,  indicative  of  the 
reciprocal  resistance,  would  be  altered. 

Conclusions 

1.  The  overall  effective  thermal  resistance  of  a  cor¬ 
rugated  flberboard  package  Is  not  constant  but  will 
vary  from  about  0.18  K»mVW  for  a  kraft  llnerboard 
package  (0.33  mm  thick)  to  about  1.18  K»mVW  for  a  cor 
rugated  flberboard  box  having  walls  51  mm  thick. 

2.  The  thermal  resistance  of  the  flberboard  Is  not  a 
reasonable  approximation  of  the  total  resistance.  The 
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resistances  of  the  box-product  intsrtscs  and  of  the 
outer  boundary  air  layer  are  significant,  adding 
0.18-0.36  K»m’)W  to  the  board  resistance.  For  normal 
box  thicknesses  the  total  resistance  will  be  3-5  times 
that  of  board  alone. 

3.  For  a  board  thickness  of  less  than  20  mm  a 
distinct  difference  in  heating  and  cooling  behavior  ex¬ 
ists  at  very  low  bulk  air  velocities  with  heating  periods 
being  significantly  greater  than  cooling  periods  on  the 
same  package. 


Appendix  I— Theory 
General 

An  energy  balance  on  a  differential  element  within  a 
solid  shows  that  heat  transfer  must  satisfy  the  basic 
differential  equation: 

dT  d*T _ k  d»T 

dt  °  dl*  pc  at1 

V  (D 


where 

k  =  thermal  conductivity  of  the  solid 
P  =  density  of  the  solid 
Cp  =  specific  heat  of  the  solid 
a  -  thermal  diffusivity  of  the  solid  ( =  k/Pcp  by 
definition) 

t  =  time  of  heat  transference 
dT  =  differential  temperature  in  the  solid 
df  z  differential  distance  at  the  point  of  differential 
temperature 
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Also,  if  there  is  a  non-negligible  thermal  resistance  be¬ 
tween  the  surface  of  the  solid  and  its  environment,  the 
solution  to  equation  (1)  must  simultaneously  satisfy  the 
boundary  equation 

q  «  UAT  (2) 

where 

q  *  heat  flux  at  the  surface 

U  *  overall  heat  transfer  coefficient 

AT  *  temperature  difference  across  the  boundary 

resistance 

Consider  a  solid  which  is  a  nonheat-produclng, 
homogeneous  rectangular  parallelepiped  at  a  uniform 
temperature,  T0.  If  it  Is  exposed  to  a  uniform 
temperature,  TA  (#  T0),  then  for  the  rectangular  coor¬ 
dinate  system  originating  at  the  center  of  ttie  solid  (fig. 
4),  the  solution*  satisfying  equations  (1)  and  (2)  has 
been  found  to  be 


*  mm.  A.  a.,  W.  £  Stawart,  and  E.  N.  Uqmfoot  Transport  Phaoomana. 
1SS0,  p.  362.  Wtlay,  N  Y 

'  Ttw  actus*  solution  is  an  inline#  aortas  solution  smtch  altar  suttlclant 
lima  will  Hava  nagHgMta  tarms  altar  tea  first  tarm  aucit  mat  may  can 
Sa  omtttad.  laavtno  aquation  (3). 


where 


TA  -  T(x,y,z) 

T .  T 

=  unaccomplished  temperature 

'A  ~  To 

change,  Tr(x,y,z) 

ta 

=  bulk  environment  temperature 

T(x,y,z) 

=  temperature  at  a  point  P(x,y,z)  in  the 

solid 

To 

ha,b,c 

=  initial  temperature  of  the  solid 
=  U/k  for  the  mutually  perpendicular 
sides  A,B,C 

=  positive  first  roots  of  the  equation 

*a,b,c 

Hl  =  Aj_  t®n  ^  “  A,B,C 

X,Y,Z 

s  axial  distances  to  the  point  P(x,y,z) 

A.B.C 

s  axial  distances  from  the  center  to 

the  surfaces  of  the  sides 

e 

s  natural  exponential 

t 

s  elapsed  exposure  time  to  TA 

Application 

A.  General  Technique  for  Finding 
Overall  Thermal  Resistance,  R, 
at  Surface  of  Solid 

Simultaneously  measure  the  outside  air  temperature 
together  with  the  solid  temperature  at  two  points  along 
the  same  axis  and  on  the  same  side  of  the  origin.  Then, 
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algebraic  manipulation*  of  equation  (3)  gives 

TA  -  T(x„0,0)  Cos  Ax,  X, 

Ta  -  T(x„0,0)  "  Cos  Ax,  X,  (4) 

at  the  time  of  the  temperature  measurement. 

Let  the  two  points  in  the  solid  be  P(A,0,0)  and 
P(0,0,0),  P(0,B,0)  and  P<0,0,0),  or  P(0,0,C)  and 
P(0,0,0).  Equation  (4)  becomes 


This  is  a  linear  relationship.  The  thermal  diffusivity  will 
be  the  negative  of  the  slope  when  fnTrfO)  is  plotted 
against  (A^  +  Ag  +  A^)t  (eq.  8a). 

This  theory  and  these  techniques  for  determining  the 
overall  heat  transfer  coefficient  and  thermal  diffusivity 
have  been  known  and  applied  for  several  decades,  but 
have  required  extensive  time  to  collect  and  analyze  the 
data.  Using  today's  electronic  data  collection  devices 
coupled  with  computer  processing,  minimal  time  and 
effort  is  needed. 


- =  Cos  A|_L 

Ta  -  T(0)  (5a) 

or  rearranging 

AL  =  Cos"'  [(TA  -  TS)/(TA  -  T(0))]/L  (5b) 

where 

A(_  =  root  of  equation  H[_  =  Al  tan  AlL 

Ta  =  bulk  air  temperature  on  side  of  the  solid 
Ts  =  surface  temperature  at  the  axis  of  the  solid 
T(O)  =  center  temperature  of  the  solid 
L  =  distance  from  center  of  the  solid  to  the  sur¬ 
face  along  the  axis 

From  the  definition  of  A  and  by  knowing  the  thermal 
conductivity  of  the  solid,  one  can  calculate  the  heat 
transfer  coefficient,  U,  and  the  thermal  resistance,  R: 

H(_  =  A|_  tan  A|_L 

U  =  kHL  =  kALtan  ALL  =  R~' 

B.  General  Technique  for  Finding 
Thermal  Diffusivity, «, 
of  Solid 

By  letting  P(x,y,z)  =  P(0,0,0),  in  equation  (3),  cos  AlL 
=  1.  Taking  the  natural  logarithm  of  equation  (3)  and 
condensing  one  gets: 

Ta  -  T(O) 

WO)  -  tn(— -)  =  C,  -  (Aa  ♦  A’g  ♦  Ab>-t8a) 
=  C,  -  (ZAJ  x  a)  t  (8b) 


C.  Modification  of  General  Techniques 
for  Determining  Resistance  and 
Diffusivity  in  a  Corrugated 
Flberboard  Package 

A  homogeneous  product  packed  tightly  in  a  corrugated 
fiberboard  box  can  be  described  as  a  solid  rectangular 
parallelepiped.  It  will  have  three  surface  thermal 
resistances  in  series:  (1)  R„  an  inner  air  resistance  be¬ 
tween  box  and  solid  arising  from  imperfect  contact;  (2) 
R„  the  resistance  due  to  the  boxboard  itself;  and  (3)  R„ 
an  outer  air  resistance  between  the  box  and  the  bulk 
air.  The  composite  resistance  is  related  to  the  overall 
heat  transfer  coefficient,  U: 


R,  +  R,  +  R,  ^overall  (9) 

To  apply  the  theory  outlined,  the  temperature  at  the  in¬ 
terface  between  the  inner  air  layer  and  the  product  sur¬ 
face  is  needed— a  measurement  almost  impossible  to 
achieve.  A  thermocouple  sandwiched  between  the  pro- 
(6)  duct  and  the  corrugated  container  would  only  give  the 
average  temperature  of  the  inner  air  layer.  Subsequent 
(T)  application  of  equations  (5b)  and  (8)  would  only  give  ap¬ 
parent  heat  transfer  coefficients  and  thermal  dif- 
fusivities  and  not  necessarily  the  true  values. 

If  InTrfO)  is  plotted  against  time  (eq.  8),  the  slope  of  the 
line  would  be  (ZAJ  x  o).  This  slope  will  be  nearly  cons¬ 
tant  since  both  time  and  inTf(0)  are  known  accurately 
and  jfnTrfO)  is  independent  of  the  interface.  Therefore,  it 
is  obvious  that 

(TA*  x  o)true  =  (TA*  x  “^apparent  0®) 

From  equation  (10),  then, 


where 


C,  =  constant  a  intercept 

Aa,8,C  -  A-values  for  mutually  perpendicular  sides 


Ta  -  T(O) 
Ta  -  T0 


unaccomplished  temperature  change  at 
center  of  solid,  Tr<0) 


'  Caralaw,  H.  S.,  and  J.  C.  Jaagr  r  conduction  of  Haat  In  Solids.  1999, 
p.  186.  Oxford  Press,  London. 


(TA*>  ,rUe 


(TA*  *  <*)  apparent 
“  true 


(11) 


Thus,  for  determining  the  true  heat  transfer  coefficient 
equation  (8b)  can  be  applied  directly  to  determine  the 
term  (ZA1  x  oiapparent  with  division  by  the  true  diffusivi¬ 
ty;  or  one  can  apply  equation  (5b)  then  equation  (8a), 
and  calculate  the  term  (ZA*  x  o)apparent  with  subse¬ 
quent  division  by  the  true  diffusivity.  While  the  former 
technique  is  straightforward,  the  latter  method  allows 
one  to  determine  the  magnitude  of  the  original  approx¬ 
imation  of  the  heat  transfer  coefficient  and  diffusivity. 
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For  this  study  the  latter  technique  was  selected  as  the 
true  diffusivity  can  be  determined  through  extrapolation 
of  the  apparent  values. 

Appendix  II— Experimental  Procedure 
System  Description 

A.  The  solid  used  was  commercial  454-gram  flat- 
packs  of  margarine  consisting  of  four  113-gram  sticks 
wrapped  in  greaseproof  paper  and  packaged  in  a  paper- 
board  box.  This  solid  is  relatively  homogeneous  and  re¬ 
tains  its  physical  properties  essentially  constant  over 
the  temperature  range  employed.  Margarine  is  also 
representative  of  actual  commercial  loads.  The  flat- 
packs  allow  for  variable  stacking  patterns.  In  this 
study,  four  packages  were  stacked  to  form  a  cube  hav¬ 
ing  approximately  130-mm  sides. 

B.  Corrugated  boxes  were  fabricated  from 
205-127C-205  corrugated  sheetstock.  Inner  flaps  were 
cut  to  meet  tightly,  and  outer  flaps  were  cut  back  to 
just  form  an  edge  seal  (flap  length  of  12.7  mm).  Box 
thickness  at  the  center  of  the  panels  was  equal  on  all 
sides.  For  greater  board  thicknesses,  larger  boxes  had 
smaller  boxes  nested  inside,  with  filler  sheets  of  cor¬ 
rugated  employed  in  the  voids. 

C.  The  environment  was  provided  by  two  adjacent 
temperature-humidity  rooms:  one  at  23°  C,  50  percent 
relative  humidity  and  the  other  at  -26°  C.  An  insulated 
window  between  the  two  rooms  allowed  access. 

D.  Air  flow  was  provided  by  a  fan  which  was  adjusted 
such  that  the  air  stream  would  strike  the  package  nor¬ 
mal  to  its  front.  Air  velocities  were  measured  at  the 
front  face  of  the  package  and  the  fan  speed  was  ad¬ 
justed  to  the  desired  level  accordingly. 

E.  Thermocouples  were  located  at  the  center  of  the 
margarine,  at  the  centers  of  the  six  interfaces  between 
margarine  and  box,  and  in  the  bulk  air  25.4  mm  from 
each  surface,  all  along  the  axes. 

F.  The  package  rested  on  12.7-mm-diameter  dowels 
set  in  a  framework  supporting  the  thermocouples  used 
to  measure  air  temperatures.  The  framework  with 
package  was  passed  between  rooms  without  disturbing 
the  relative  positions  of  the  thermocouples. 

G.  A  programmable  data  acquisition  system  would 
periodically  scan  the  thermocouples.  Combined  scan- 
output  time  was  approximately  20  seconds  with  actual 
scanning  time  being  much  less. 


Procedure 

A.  Data  Collection 

Thermal  equilibrium  was  established  in  the  package. 

An  initial  scan  just  prior  to  starting  the  run  gave  the  in¬ 
itial  temperature,  T0.  With  the  fan  on,  the  package  and 
framework  were  quickly  subjected  to  a  step  increase  in 


bulk  air  temperature  by  being  transferred  to  the  adja¬ 
cent  room  through  the  access  panel  in  the  common 
wall.  The  data  acquisition  system  was  activated  at  time 
zero.  Temperature  scans  were  taken  until  the  center 
temperature  approximated  the  bulk  air  temperature. 
Scan  intervals  were  adjusted  periodically  to  reflect  rate 
of  temperature  change,  and  approximately  100  scans 
per  test  were  collected.  When  the  package  temperature 
had  again  equilibrated,  the  frame  and  package  were 
brought  back  into  the  original  room  and  a  second  set 
of  data  was  recorded.  By  this  method  both  heating  and 
cooling  data  were  collected  for  each  test  package  and 
air  velocity. 

B.  Data  Processing 

Initial  data  processing  consisted  of  calculating 

TA  ~  TS 

Ta  -  T(O)  and  inTr(0)  for  use  in  equations  (5b,  6,  and 
8).  These  values  were  determined  for  each  side  of  the 
package  and  for  each  time  scan. 

Since  the  study  was  intended  to  determine  an  average 
value  for  U  which  would  be  indicative  of  box  behavior, 
the  next  processing  step  integrated  the  time-smoothed 
values  of  An  and  calculated  for  each  side  of  the 
package  an  average  An  according  to  the  formula  (fig. 
3A,B): 

-  *o  An(t)At 

An  =  - f - — 

to  ^  (12) 

From  An,  Hn  wasjletermined  by  equation  (6).  The 
overall  apparent  J^and  H  for  the  box  was  calculated  by 
averaging  the  six  An  and  Hn  values. 

Calculation  of  the  apparent  thermal  diffusivity  was 
more  complicated.  In  order  to  determine  the  slope 
(eq.  8),  the  value  of  IAJ  =  (A*  +  Ay  +  A|)  needed  to  be 
known;  but  for  each  side  x  and  its  corresponding  Ax 
there  are  four  combinations  of  Ay  and  Az.  Therefore  all 
four  IA-  were  determined  and  averaged.  Once  this  was 
done  fnTr(0)  could  be  plotted  verses  (XA-)t.  An  algorithm 
was  then  applied  to  find  the  linear  region  of  the  data. 
From  here  the  slope  was  found.  This  process  also  was 
repeated  for  each  side  (fig.  5).  The  six  values  of  ap¬ 
parent  thermal  diffusivity  were  averaged  as  the  solid 
was  not  perfectly  homogeneous. 

The  final  output  gave  the  apparent Jhermal  diffusivity, 
and  the  apparent  values  of  An  and  Hn  for  each  of  the 
six  sides  of  the  package  followed  by  the  overall 
averages  of  the  six  sides,  including  standard  deviations 
and  coefficients  of  variation  among  the  sides.  These 
values  constituted  the  data  used  for  the  initial  analysis 
of  the  effects  due  to  air  velocity  and  board  thickness. 

C.  Determination  of  Thermal 
Conductivity  of  Margarine 
Flatpacks 

The  data  processing  yields  values  of  H,  which  is  the  ap- 
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Figure  6  —  Thermal  conductivity  ot  Poxed  margarine  using  a  guarded  hot 
plate  apparatus  Standard  deviation  eguaieo  0  005  W'Mm  Coet 
ticient  of  variation  was  2  8  percent 
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parent  overall  heat  transfer  coefficient  divided  by  the 
thermal  conductivity  of  the  margarine,  U/k.  To  deter¬ 
mine  U,  it  was  necessary  to  first  determine  k.  This  was 
accomplished  using  a  commercial  thermal  conductivity 
tester  calibrated  against  a  National  Bureau  of  Stan¬ 
dards  glass  fiber  standard  and  based  on  the  ASTM 
C-518  procedure.  Flatpacks  of  margarine,  31.8  mm 
thick,  were  used  as  the  test  sample.  Two  temperature 
differentials  were  employed:  9.5°  C  and  19.2°  C.  The 
data  are  presented  in  figure  6.  The  scatter  was  about 
equal  for  both  temperature  gradients.  For  the 
temperature  range  in  this  study  the  thermal  conductivi¬ 
ty  only  varied  about  2.5  percent  from  0.1744  W/K»m  at 
-26°  C  to  0.1786  W/K*m  at  27°  C.  This  permitted  the 
use  of  the  average  conductivity  which  was  0.176  W/K*m 
with  a  standard  deviation  of  0.005  W/K»m  and  a  coeffi¬ 
cient  of  variation  of  2.81  percent. 

D.  Determination  of  True  Thermal 
Diffusivity  of  Margarine 

To  obtain  the  true  thermal  diffusivity  the  data  for 
margarine  alone  (no  package)  was  employed  (fig.  7). 
Only  the  cooling  data  was  used,  as  the  thermal 
resistance  was  lower  than  for  heating.  Extrapolation  to 
zero  reciprocal  air  velocity  gives  the  diffusivity  at  in¬ 
finite  air  velocity  where  the  outer  boundary  resistance 
would  theoretically  cease  to  exist.  A  linear  regression 
of  the  three  points  gives  an  intercept  value  of  9.406  x 
10' 8  mJ/s.  This  value  was  assumed  to  be  the  true  ther¬ 
mal  diffusivity  of  the  product.  Applying  equation  (11)  in 
the  following  form  gives  the  true  A-value: 


+  **  +  X‘z)  *  °1  parent 

(*x  +  Ay  +  *z)  -  — - - 

°  true 


(13) 


Since  the  composite  A-value  was  calculated  to  represent 
the  package  as  a  whole. 
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figure  7  —  Determination  ot  true  thermal  diffusivity  ot  packaged  margarine 
Extrapolation  ot  data  to  0  <misi  ’  gives  thermal  diffusivity  for 
a  perfectly  mixed  heat  source  or  sink  which  is  the  true 
thermal  diffusivity  of  the  product 

(M  149721-1) 


Ax  =  Ay  =  Az  =  A,  the  true  average  A-value 
and  j [^  =  Ay  =  Az  =  A,  the  apparent  average  A-value 
and  equation  (13)  becomes 


3  Mrue  = 


(3  A2  x  o)  apparent 


“true 


or, 


_  Tv  <  °  apparent 
Afrue  -  *  x  t  I 

“true 


(14a) 


(14b) 
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E.  Final  Data  Calculation 

Using  equations  (14b  and  6),  plus  the  output  from  the 
initial  data  processing,  plus  the  true  diffusivity,  one 
gets  the  desired  resistance  for  the  various  board 
thicknesses  and  air  velocities,  The  resistance  values  at 
each  thickness  were  then  plotted  versus  reciprocal  air 
velocity  and  linear  regression  lines  were  calculated  us¬ 
ing  0.09  m/s  as  an  estimate  of  0  m/s  (table  2).  These 
data  were  extrapolated  to  0  (m/s)  '  for  the  infinite  air 
velocity  values  for  table  1  and  figure  1. 


Appendix  III— Hypothesis  Concerning  Differences 
In  Heating  and  Cooling  Behavior 

Situation  1— Heating  of  Solid  Product 

When  heating  the  product,  thermal  energy  moves  from 
the  bulk  air  through  a  semistagnant  air  film  layered  be¬ 
tween  the  bulk  air  mass  and  the  box.  This  boundary 
layer  will  have  a  finite  effective  thickness  and  a  ther¬ 
mal  gradient,  and  thus  also  a  density  gradient.  As  the 
bulk  air  velocity  increases,  this  film  thickness 
decreases  and  correspondingly  its  thermal  resistance 
decreases.  This  is  evident  on  figure  1  where  line  II  is 
less  than  line  III.  Also,  as  the  board  thickness  is  re¬ 
duced.  the  overall  thermal  resistance  decreases  and 
the  overall  thermal  flux  increases.  In  the  absence  of 
sufficient  air  mass  movement  to  re-supply  the  energy 
draining  from  the  bulk  air.  that  part  of  the  bulk  air  adja¬ 
cent  to  the  outer  edge  of  the  stagnant  air  layer  will 
cool  and  itself  become  part  of  the  stagnant  layer,  i.e.. 
the  stagnant  layer  thickens  and  its  resistance  value  in¬ 
creases.  This  is  evident  in  figure  1.  line  lll-a.  where  the 
difference  in  resistance  between  lines  lll-a  and  II  in¬ 
creases  even  though  the  overall  resistance  (line  lll-a)  is 
decreasing. 

On  the  other  hand,  an  increased  air  velocity  maintains 
the  necessary  bulk  thermal  energy  source.  The  in¬ 
creased  velocity  will  stir  the  outer  edge  of  the  stagnant 
air  layer  making  it  thermally  more  uniform  with  the  bulk 
air.  thereby  effecting  a  reduction  in  the  thickness  of  the 
remaining  air  layer  and  its  corresponding  resistance. 
The  overall  resistance,  line  III.  would  converge  on  line 
II.  Examination  of  the  data  of  table  1  for  any  board 
thickness  and  for  the  different  air  velocities  shows  this 
effect.  The  key  point  is  that  a  change  in  thickness  of 
the  outer  stagnant  air  boundary  layer  due  to  the 
presence  or  nonpresence  of  air  flow  has  a  significant 
effect  on  overall  resistance. 


Situation  2— Cooling  the  Solid  Product 

For  the  opposite  case,  cooling,  an  entirely  different 
mechanism  may  be  taking  place.  With  the  outward  heat 
flux  there  will  be  a  reversal  of  the  temperature  gradient 
across  each  resistance  zone  relative  to  heating  and  a 
corresponding  reversal  of  the  density  gradient  of  the 
outer  boundary  air  layer.  This  situation  places  the  lower 
density  air  adjacent  to  the  board.  As  in  the  case  of 
heating,  at  the  greater  board  thicknesses  the  thermal 


flux  is  sufficiently  low  that  the  air  film  will  remain 
relatively  undisturbed  and  the  resistance  behaves  as 
for  heating.  For  the  lower  range  of  board  thickness  the 
overall  resistance  is  lower,  as  it  was  for  heating 
However,  the  higher  outward  thermal  flux  accompany 
ing  this  lower  resistance  may  cause  thermal  hot  spots 
to  develop  at  high  points  of  the  board-air  interface 
Because  of  the  heterogeneous  nature  of  the  surface, 
these  hot  spots,  located  at  the  fiber  tips,  may  be  the 
source  of  thermal  convection  currents.  These  currents 
may  tend  to  disrupt  the  air  film  bond  at  the  board  sur 
face,  thereby  allowing  the  layer  to  be  virtually  stripped 
away  by  even  a  minimal  bulk  air  flow,  decreasing  the 
overall  resistance  As  the  board  thickness  is  reduced 
further  the  rate  of  convective  transfer  will  increase  due 
to  increased  flux  and  the  resistance  III  b  of  the  stag 
nant  layer  decreases  relative  to  resistance  III  a  If  the 
flux  is  sufficiently  high,  the  entire  stagnant  layer  nearly 
ceases  to  exist  and  the  resistance  values  approach 
those  for  infinite  air  velocity  line  II 

To  show  that  there  is  actually  a  mechanistic  difference 
between  heating  and  cooling,  one  need  only  examine 
the  time  histories  (fig.  3)  for  a  no-air-flow  single-wall 
box.  The  value  "A"  is  reflective  of  the  overall  heat 
transfer  coefficient  and  thus  the  reciprocal  of  the 
overall  resistance.  At  the  beginning  of  the  heating  run 
(fig.  3-  B)  there  is  a  large  inward  thermal  flux  due  to  the 
large  overall  temperature  difference.  The  high  flux 
disrupts  the  air  layer  almost  completely,  giving  high 
heat  transfer  efficiency,  shown  by  the  high  values  of 
"k"  at  the  beginning  of  the  run  As  the  overall 
temperature  difference  between  bulk  air  and  product 
decreases  with  time,  so  does  the  overall  flux.  The  dis¬ 
turbance  of  the  outer  air  layer  lessens  and  the  layer  ef¬ 
fectively  thickens,  thereby  increasing  its  resistance  and 
decreasing  the  flux  even  further.  This  is  evidenced  by 
the  rapid  decrease  in  "A."  For  cooling  (fig.  3-A)  the  layer 
is  sufficiently  thin  due  to  the  postulated  convective  ef 
feet  that  the  resistance  stays  neariy  constant  and  at  a 
low  level,  as  is  evidenced  by  the  higher  value  of  "A" 
throughout  the  run.  The  shapes  of  these  curves  are 
typical  of  all  runs  made,  varying  only  in  accordance 
with  the  mechanism  proposed. 
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While  the  mechanical  properties  of  corrugated 
fiberboard  are  known,  little  knowledge  has  been 
available  on  its  thermal  properties.  Here  the  thermal 
resistance  of  a  container  was  determined  as  a  function 
«f  air  velocity  and  board  thickness,  and  an  anpropriate 
design  curve  is  presented. 
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